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ABSTRACT 



Aims. The aim of the present paper is to identify and study the properties and galactic content of groups and clusters in the GOODS- 
South field up to z ~ 2.5, and to analyse the physical properties of galaxies as a continuous function of environmental density up to 
high redshift. 

Methods. We use the deep (z 850 ~ 26), multi-wavelength GOODS-MUSIC catalogue, which has a 15% of spectroscopic redshifts 
and accurate photometric redshifts for the remaining fraction. On these data, we apply a (2+l)D algorithm, previously developed by 
our group, that provides an adaptive estimate of the 3D density field. We support our analysis with simulations to evaluate the purity 
and the completeness of the cluster catalogue produced by our algorithm. 

Results. We find several high density peaks embedded in larger structures in the redshift range 0.4-2.5. From the analysis of their 
physical properties (mass profile, M 2 m, o" v , L x , if - B vs. B diagram), we derive that most of them are groups of galaxies, while two 
are poor clusters with masses of few times 10 14 M G . For these two clusters we find, from the Chandra 2Ms data, an X-ray emission 
significantly lower than expected from their optical properties, suggesting that the two clusters are either not virialised or gas poor. 
We find that the slope of the colour magnitude relation, for these groups and clusters, is constant at least up to z ~ 1. We also analyse 
the dependance on environment of galaxy colours, luminosities, stellar masses, ages and star formations. We find that galaxies in high 
density regions are, on average, more luminous and massive than field galaxies up to z ~ 2. The fraction of red galaxies increases with 
luminosity and with density up to z ~ 1.2. At higher z this dependance on density disappears. The variation of galaxy properties as a 
function of redshift and density suggests that a significant change occurs at z ~ 1.5-2. 

Key words. Galaxies:distances and redshift - Galaxies: evolution - Galaxies: high redshift - Galaxies: clusters: general - Galaxies: 
fundamental parameters - (Cosmology:) large-scale structure of Universe 



1. Introduction 

The study of galaxy clusters and of the variation of galaxy 
properties as a function of the environment are fundamental 
tools to understand the formation and evolution of the large 
scale structures and of the different galaxy populations, observed 
both in the local and in the high redshift Universe. The ef- 
fects of the environment on galaxy evolution have been stud- 
ied at progressiv ely higher redshifts through the analysis of sin- 
gle clusters (e.g [Treu et alj|2003t iNakata et alfeOOl iTran et alJ 



12001 iMei etafl boOot iMenci et all l2Q08h iRettura et alJ I2008I) 
as well as studying the variation of galaxy colours, mor- 
phologies and other physical parameter s as a funct ion of 
projected or 3-dimensional density (e.g. iDressler et alJ 



y i 

Blanton et all |2005t ICucciati et atl l2006l: ~ ICooper et all 



1997 



2007 



Elba z et al. [ 120071) . Moreover, the analysis of cluster properties 



at different wavelenghts provides interesting insights into the 
matter content and evolutionary histories of these structure s 
dLubin et alj|2004t iRasmussen et alj|2006t iPopesso et aT]|2007l) . 
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A variety of survey techniques have proved effective at find- 
ing galaxy clusters up to z~ 1 and beyond. X-ray selected 
samples at z > 1 pr obe the most massive and dynamically re- 
laxed systems (e.g . . iMaughan et al.ll2Q04t IStanford etai1 [2006; 
iBremeret alj |2006: LidmanetaL|[2008). Large- area multicolour 
surve ys, such as the red-sequence survey (e.g. iGladders & Yed 
2005), have collected samples of systems in a range of evo- 
lutionary stages. The mid-IR cameras on board of the Spitzer 
and Akari satellites has extended the range and power of multi- 
colour s urveys, producing confirmed and candidate clusters up to 
z ~ 1.7 (IStanford et alj|2005t lEisenhardt et alj|2008t iGoto et ail 
2008). However most of the previous techniques present some 
difficulties in the range 1.5 < z < 2.5, where we expect to ob- 
serve the formation of the red sequence and the first hints o f 
colour segregation (ICucciati et al.1 120061: iKodama et al.l 12007). 
Searching for extended X-ray sources becomes progressively 
more difficult at large distances, because the surface bright- 
ness of the X-ray emission fades as (1 + z) 4 . The sensitiv- 
ity of surveys exploiting the Sunyaev-Zeldovich (SZ) effect is, 
at pre sent, not sufficient to detect any of known c lusters at 
z > 1 (ICarlstrom et al.ll2002l: Istaniszewski et al]|2008l) . Finally, 
the detection of galaxy overdensities on surveys using two- 
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dimensional algorithms, requires additional a priori assumptions 
on either galaxy luminosity function (LF), as in the Matched 
Filter algorithm (Postman et al. 1996), or relies on the presence 
of a red sequence (Gla dders & Yedl2000l) . Biases produced by 
these assumptions can hardly be evaluated at high redshift. 

In this context, photometric redshifts obtained from deep 
multi-band surveys for large samples of galaxies, though hav- 
ing a relatively low accuracy if compared to spectroscopic red- 
shifts, can be exploited to detect and stud y distant structures . 
In the past few years, several authors (e . g. iBotzler et aL 2004; 



van Breukelen et all l2006| IScoville et al.l 120071: [Z atloukal et alJ 
20071: iMazure et alj|2007t lEisenhardt etaTl 2008) have devel- 



oped or extended known algorithms to take into account the 
greater redshift uncertainties. We have developed a new algo- 
rithm, that uses an ada ptive estimate of the 3D density field, as 
described in detail in iTrevese et alJ J2007). This method com- 
bines galaxy angular positions and precise photometric redshifts 
to estimate the galaxy number-density and to detect galaxy over- 
densities in three dimensions also at z > 1, as described in Sect. 
3. 

Our first application to the K20 survey (Cim atti et al.ll2 002) 
detected two clusters at z~0.7 and z~l (ITrevese et alJ 12007). 
previously identifi ed through spectroscopy dGilli et al.l [2003; 
Adami et al. 2005). We then appli ed the algorithm to the much 
larger GOODS-South field, and in lCastellano et al.l (120071) . else- 
where C07, we reported our initial results, i.e. the discovery of a 
forming cluster of galaxies at z~ 1.6. 

In this paper we present the application of the algorithm 
to the entire GOODS-South area (~ 143 arcmi n 2 ), using the 
GOODS-MUSIC catalogue (iGrazian et al.ll2006al) up to z ~ 2.5, 
to give a comprehensive description of the large scale structures 
in this field, with a detailed analysis of the physical properties of 
each high density peak. We also study the physical properties of 
galaxies as a function of environmental dens ity up to redshift 2.5 , 
higher than previo us similar studies (e.g.. ICucciati et ail [2006; 
ICooper et alj|2007l) . 

To validate our technique, we analysed the completeness and 
purity of our cluster detection algorithm, up to z ~ 2.5, through 
its application to a set of numerically simulated galaxy cata- 
logues. Besides allowing an assessment of the physical reality 
of the structures found in the GOODS field, this analysis pro- 
vides the starting point to test the reliability of the algorithm in 
view of our plan to apply it to photometric surveys of similar 
depth but covering much larger areas. 

The paper is organised as follows: in Sect. 2, we describe the 
basic features of our dataset. In Sect. 3, we summarise the basic 
features of the (2+l)D algorithm used in our analysis, and com- 
pare it with other methods based on photometric redshifts. In 
Sect. 4, we show the results of the application of our method to 
simulated data. In Sect. 5, we present the catalogue of the struc- 
tures detected and the derived physical properties. In Sect. 6, we 
study the colour magnitude diagrams of the detected structures. 
In Sect. 7, we analyse the physical properties of galaxies as a 
continuous function of environmental density. 

All the magnitudes used in the present paper are in the AB 
system, if not otherwise declared. We adopt a cosmology with 
Q A = 0.7, Q M = 0.3, and H = 70 km s 1 Mpc" 1 . 



2. The GOODS-MUSIC catalogue 

We used the multicolour GOODS-MUSI C catalogue (GOOD S 
Multicolour Southern Infrared Catalogue: lGrazian et al . 2006a). 
This catalogue comprises information in 14 bands (from U band 



to 8//m) over an area of about 143.2 arcmin 2 . We used the zsso- 
selected sample (z%so ~ 26), that contains 9862 galaxies (after 
excluding AGNs and galactic stars). About 15% of the galax- 
ies in the sample have spectroscopic redshift, and for the other 
galaxies we used photometric redshifts obtained from a standard 
y 2 minimisation oy er a large set of spectral models (see e.g., 
Font ana et al.l l2000). The accuracy of the photometric redshift is 
very good, with a r.m.s. of 0.03 for the Az/(1 + z) distribution up 
to reds hift z — 2. For a detail ed description of the catalogue we 
refer to lGrazian et al.l (l2006al) . 

The method we applied to estimate the rest-frame magni- 
tudes and the other physical par ameters (M, SFR, a ge) is de- 
scribed in previous papers (e.g., iFontana et al.l 2006). Briefly, 
we use a x 2 minimisation analysis, comparing the observed 
SED of each galaxy to synthetic templates, and the redshift 
is fixed during the fitting pr ocess to the spectrosc opic or pho- 
tometric redshift derived in iGrazian et al.l d2006al) . The set of 
templates is computed w ith standard spectral synthesis models 
(iBruzual & Charlotll2003l) . chosen to broadly encompass the va- 
riety of star formation histories, metallicities and extinctions of 
real galaxies. For each model of this grid, we compute the ex- 
pected magnitudes in our filter set, and find the best-fitting tem- 
plate. From the best-fitting template we obtain, for each galaxy, 
the physical parameters that we use in the analysis. Clearly, the 
physical properties are subject to uncertainties and biases related 
to the synthetic libraries used to fit the galaxy SEDs; however, 
as shown in lFontana et al.l (|2006), the extension of the SEDs to 
mid-IR wavelengths with IRAC tends to reduce the uncertain- 
ties on the derived stellar masses. For a detailed analysis of the 
uncerta inties on the physical properties we refer to our previous 
papers (IFontana et alJl2006t IGrazian et al.ll2007l) . In the present 
work we also make use of the 2Ms X-ra y observation o f the 
Chandra Deep Field South presented by iLuo et all (120081) and 
of the catalogu e of VLA radio sources (1.4 GHz) on the CDFS 
compiled bv lMiller etail (120081) . 

3. The (2+1 )D algorithm for the density estimation 

To estimate a three dimensional density, we developed a method 
that combines the angular position with the photometric red- 
shift of each objec t. The algorithm is described in detail in 
ITrevese et al. (2007): here we outline its main features and, in 
the next section, we present the simulations used to estimate its 
reliability. 

The procedure is designed to automatically take into account 
the probability that a galaxy in our survey is physically associ- 
ated to a given overdensity. This is obtained by computing the 
galaxy densities in volumes whose shape is proportional to po- 
sitional uncertainties in each dimension (a, 6 and z). 

First, we divide the volume of the survey in cells whose ex- 
tension in different directions (Aa, A6, Az) depends on the rel- 
evant positional accuracy and thus are elongated in the radial 
direction. We choose the cell sizes small enough to keep an ac- 
ceptable spatial resolution, while avoiding a useless increase of 
the computing time. We adopt Az = 0.025 (radial direction) and 
Aa = AS ~ 2.4 arcsec in transverse direction, the latter value 
corresponds to ~ 30, 40 and 60kpc (comoving), respectively at 
z ~ 0.7, ~ 1.0 and ~ 2.0. 

For each cell in space we then count neighbouring ob- 
jects in volumes that are progressively increased in each di- 
rection by steps of one cell, thus keeping the simmetry im- 
posed by the different intrinsic resolution. When a number n 
of objects is reached we assign to the cell a comoving density 
p = n/V„, where V n is the comoving volume which includes 
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the n-nearest neighbours. Clusters would be better characterized 
by their proper density since they already decoupled from the 
Hubble flow, however we notice that the average uncertainty on 
photometric redshifts, that grows with redshift as (1 + z), forces 
us to measure densities in volumes that are orders of magnitude 
larger than the real volume of a cluster, even at low-z. Thus we 
decided to measure comoving densities, that have the further ad- 
vantage of giving a redshift-independent density scale for the 
background. We fix n — 15 as a trade off between spatial reso- 
lution and signal-to-noise ratio. Indeed, through simulation de- 
scribed in Sect. |4] we verified that a lower n would greatly raise 
the high frequency noise in the density maps, thus increasing 
the contamination from false detections in the cluster sample, 
even at low redshift ('purity' parameter in Tab. [TJ. In the den- 
sity estimation, we assign a weight w(z) to each detected galaxy 
at redshift z, to take into account the increase of limiting ab- 
solute magnitude with increasing redshift for a given apparent 
magnitude limit. We choose w(z) - l/s(z), where s(z) is the 
fraction of objects detected with respect to a reference redshift 
z c below which we detect all objects brighter than the relevant 
M c = M lim {ZcY 



s{z) = 



f M ^®(M)dM ' 



(1) 



where <D(M) is the redshift dependent galaxy luminosity 
function computed on the same GOODS-MUSIC catalogue 
dSalimbeni et al.ll2008l) . M/ !m (z) is the absolute magnitude limit 
at the given redshift z, corresponding to the apparent magni- 
tude limit m;,m of the survey, which depends on the position (see 
iGrazian et al. 2006a). We use this correction to obtain a density 
scale independent of redshift, at least to a first approximation. In 
computing M;, m (z), we use K - and evolutionary - corrections for 
each object computed with the same best fit SEDs used to de- 
rive the stellar masses, the rest-frame magnitudes, and the other 
physical properties. 

We apply this algorithm to data from the GOODS-MUSIC 
catalogue, in a redshift range from z ~ 0.4 to z ~ 2.5, where 
we have sufficient statistic. We perform this analysis selecting 
galaxies brighter than Mb = -18 up to redshift 1.8 and brighter 
than Mb = -19 at higher redshift, to minimise the completeness 
correction described above, keeping the average weight w (z) be- 
low 1 .6 in all cases. 

Using this comoving density estimate we analyse the field 
in two complementary ways. First, we detect and study galaxy 
overdensities, i. e. clusters or groups (see Sect. defined as 
connected 3 -dimensional regions with density exceeding a fixed 
threshold and a minimum number of members chosen accord- 
ing to the results of the simulations (Sect. 0). In particular, we 
isolate the structures as the regions having p > p + 4<x on our 
density maps and at least 5 members. We then consider as part 
of each structure the spatially connected region (in RA, DEC, 
and redshift) around each peak, with an environmental density 
of > 2<x above the average and at least 15 member galaxies. To 
avoid spurious connections between different structures at the 
same redshift, we consider regions within an Abell radius from 
the peak. The galaxies located in this region are associated with 
each structure. We then study the variation of galaxy properties 
as a function of environmental density (Sect. |7), associating to 
each galaxy in the sample the comoving density at its position. 



3.1. Comparison with other methods based on photometric 
redshifts 

As mentioned in the introduction, other methods, based on pho- 
tometric redshifts, have been developed for the detection of cos- 
mic structures. Here we present the main differences between 
our algorithm and those which appeared most recently in the 
literature. However, a more detailed comparison is beyond the 
scope of the present work, since it would require extensive sim- 
ulations and/or the application of the different methods to the 
same datasets. 

A similar three-dim ensional approach has been proposed by 
IZatloukal et alJ (I2007I) . They select cluster candidates detecting 
excess density in the 3D galaxy distribution reconstructed from 
the photometric redshift probability distributions p(z) of each 
object. However, at variance with our method they do not adopt 
any redshift dependent correction for their estimated density, 
since they analyze only a small redshift range. As we outlined in 
the previous section, such correction is needed to provide a red- 
shift independent density scale in a deep sample as the GOODS- 

M USIC. 

iBotzler et alJ (|2004|) expanded the well kn own Friends of 
Friends (FoF) algorithm dHuchra & Gellerll982l) . to take into ac- 
count photometric redshift uncertainties. This method links to- 
gether groups of individual galaxies if their redshift difference 
and angular distances are below fixed thresholds. These thresh- 
olds depend on the photometric redshift uncertainties, which are 
greater than the average physical distance between galaxies and 
also greater than the velocity dispersion of rich clusters. This 
could induce the problem of structures percolating through ex- 
cessively large volumes. They dealt with this issue dividing the 
catalogue in redshift slices. Instead of comparing the distance 
between galaxy pairs, as done in a FoF approach, we use the 
statistical information of how many galaxies are in the neigh- 
bourhood of a given point to estimate a physical density. This 
approach can avoid more effectively the percolation problem, 
since it identifies structures from the Act density peaks whose 
extension is limited by the fixed threshold in density. 

Scoville et all (l2007l), iMazure et all 



Several authors, e. 



(12007 ). Eisenh ardt et~aTTd2008l) and Ivan Breukelen et all (2006) 



estimated the surface density in redshift slices, each with dif- 
ferent methods: the first two u se adaptive smoothing of galaxy 
counts. lEisenhardt et al.l d2008l) analyse a density map convolved 
with a wavelet kerne l, while the last auth or adopts FoF and 
Voronoi tessellation (Marin oni et alJ [2002). At variance with 
these, we prefer to adopt an adaptive 3D density estimate to con- 
sider, automatically, distances in all directions and the relevant 
positional accuracies at the same time. This approach requires 
longer computational times, but allows for an increased resolu- 
tion in high density regions where the chosen number of objects 
is found in a smaller volume with respect to field and void re- 
gions. As a consequence it also avoids all peculiar "border" ef- 
fects given by the limits of the redshift slices, and there is also 
no need to adopt additional criteria to decide whether an over- 
density, present in two contiguous 2D density maps in similar 
angular positio ns, represents the same group or not (as done for 
example by IMazure et al.ll2007l) . This clearly also depends on 
the ability of the algorithm to separate aligned structures (for a 
more detailed discussion of this see Sect. |4). 

Finally, another important difference with respect to pre- 
vious methods is in the way we use the photometric 
redshift: some a uthors used b e st fit values of photomet- 
ric redshift, e.g. Mazure et al. d2007l). while IScoville et all 
( 2007b, Ivan Breukelen et alJ d2006l) . IZatloukal et al.l d2007l) and 
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Eisenhard t et al.l ((2008) consider the full probability distribution 
function (PDF ) to take into accou nt redshift uncertainties. As 
discussed by IScoville et alj J2007), this last method could tend 
to preferentially detect structures formed by early type galaxies, 
since they have smaller photometric redshift uncertainty, thanks 
to their stronger Balmer break, when this feature is well sampled 
in the observed bands. We are less biased in this respect, since 
we consider the photometric redshift uncertainty in a conserva- 
tive way, choosing only the maximum redshift range where we 
count neighbour galaxies to associate with each cell. We took 

this range as ±2 • o> around the redshift of each cell, where 

r ii 1 

cr z = 0.03 • (1 + z) (IGrazian et al.ll2006al) is the average accu- 
racy of the photometric redshift in the range we analyze. 

4. Simulations 

We estimate the reliability of our cluster detection algorithm 
by testing it on a series of mock catalogues, designed to repro- 
duce the characteristics of the GOODS survey. These mock cat- 
alogues are composed by a given number of groups and clusters 
superimposed on a random (poissonian) field. While this is a 
rather simplistic representation of a survey, it allows us to eval- 
uate some basic features of our algorithm, without the use of N- 
bo dy simulations . We ex pand the previous simulations presented 
in iTrevese et al.l d2007l) . using a larger number of mock cata- 
logues and adopting a more consistent treatment of the survey 
completeness. For each redshift, we calculate the limiting abso- 
lute B magnitude for the two populations of "red" and "blue" 
g alaxies, defined from th e minima in the U-V vs. B distribution 
in lSalimbeni et al .1 (2008), using the average type-dependant K- 
and evolutionary corrections calculated from the best-fit SED of 
the objects in the real catalogue. We then generate an "observed" 
mock catalogue of field galaxies randomly distributed over an 
area equal to that of the GOODS-South survey. At each red- 
shift, the number of objects in the catalogue is obtained from the 
integral of the rest frame B band luminosity function <&(Mb,z) 
derived in Salimb eni et al.l d2008l) . up to the limiting absolute 
M B (z) magnitude computed as described above. 



Table 1. Completeness and Purity 



Table 2. Average distances of detected peaks from real centres. 



Mass 



Purity Un. Pairs Compl. Double Id. 









0.4 <z 


< 1.2 




M> 1 x 10 13 


M s 


100% 


6.5% 


86.2% 


4.2% 


M>2x 10 13 


M 


100% 


3.1% 


89.7% 


0% 


M> 3 x 10 13 


M e 


93.8% 


0% 


93.8% 


0%' 








1.2 <z 


< 1.8 




M> 1 x 10 13 


M e 


95.4% 


4.6% 


76.4% 


5.7% 


M> 2 x 10 13 


M s 


95.3% 


4.6% 


84.3% 


0% 


M> 3 x 10 13 


M e 


100% 


0% 


82.4% 


2.9% 


1.8 <z<2.5 


M> 1 x 10 13 


M G 


76% 


3% 


33% 


0% 


M> 2 x 10 13 


M e 


78% 


0% 


39% 


0% 


M> 3 x 10 13 


M s 


75% 


0% 


37% 


0% 



Finally, we create different mock catalogues superimposing 
a number of structures on the random fields. Given the rela- 



Redshift Interval Ar (Mpc) 



A; 



0.4 < z < 1.2 
1.2 <z< 1.8 
1.8 <z<2.5 



0.13 ±0.09 
0.13 ±0.07 
0.20 ±0.11 



0.016 ± 0.013 
0.028 ±0.019 
0.044 ±0.033 



Table 3. Separation threshold for aligned groups. 



Distance 



Separation threshold 



Redshift Projected 



1 



la. 



Ar=1.0Mpc 
Ar=1.5 Mpc 
Ar=2.0 Mpc 



> 8cr 
6a 

4(7 



> 8cr 
5a 
5<r 



2a- 



Ar=1.0 Mpc 
Ar=1.5 Mpc 
Ar=2.0 Mpc 



> 8cr 
4a 
4a 



> 80- 
5a 
5a 



3a. 



Ar=1.0 Mpc 
Ar=1.5 Mpc 
Ar=2.0 Mpc 



5a 
4a 
4a 



5a 
5a 
5a 



tively small comoving volume sampled by the survey, we ex- 
pect to find only groups and small clusters with a total mass 
M ~ 10 13 - 10 14 M and a number of members corres ponding to 
the lowest Abell richness classes dGirardi et aljl9 98a). To check 
that the performance of the algorithm does not change appre- 
ciably with a varying number of real overdensities of this kind, 
we perform three different subset of simulations. Each subset is 
based on the analysis of 10 mock catalogues, with a number of 
groups equal to the number of M > 10 13 M Q , M > 2 x 10 I3 M o 
and M > 3 x 1O I3 M DM haloes, obtained by integ rating the 
Press & Schechter function (iPress & Schechter] 1 19741) over the 
comoving volume sampled by the survey. Their positions in real 
space are chosen randomly. Cluster galaxies foll ow a King-like 
spatial distribution n(r) oc [1 + (r/r c ) 2 ]~ 3/2 (see ISarazin| [l988) 
with a typical core radius r c = 0.25Mpc. 

To take into account the uncertainty on photometric red- 
shifts, to each cluster we assign galaxy a random redshift ex- 
tracted from a gaussian distribution centred on the cluster red- 
shift Zd and having a dispersion o~ z = 0.03 • (1 + z c d- We neglect 
the cluster real velocity dispersion, which is much smaller than 
the Zphot uncertainty. We analyse the simulations in the same way 
as the real catalogue, i.e. calculating galaxy volume density con- 
sidering objects with Mb < — 18 at z < 1.8 and objects with 
M B < -19 at z> 1.8. 

We evaluate the completeness of the sample of detected clus- 
ters (fraction of real clusters detected) and its purity (fraction of 
detected structures corresponding to real ones) at different red- 
shifts (see Table [TJ. We present also the number of unresolved 
pairs (a detected structure corresponding to two real ones) and 
the number of double identifications (a unique real structure sep- 
arated in two detected ones). 

Our aim is to study the properties of individual structures and 
not, for example, to perform group number counts for cosmolog- 
ical purposes. For this reason, we prefer to choose conservative 
selection criteria in order to maximise the purity of our sample, 
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while still keeping the completeness high. We isolate the struc- 
tures as described in Sect. [3] and we consider as significant only 
those overdensities with at least 5 members in the 4cr region and 
15 members in the 2<x region. 

A structure in the input catalogue is identified if its center is 
within Ar = 0.5Mpc projected distance, and within Az = 0.1, 
from the center of a detected structure, for the low redshift sam- 
ple and Ar = O.SMpc, Az = 0.2 at high z (to account for the 
increased uncertainties in redshift and position). The results are 
reported in Table Q] We can see that the chosen thresholds and 
selection criteria allow for a high purity (~ 100%) at z < 1 .8, still 
detecting about the 80% of the real structures. At z > 1.8, given 
the greatly reduced fraction of observed galaxies, the noise is 
higher and these criteria turn out to be very conservative (there- 
fore the completeness is low) but are necessary to keep a low 
number of false detections (purity ~ 75 - 80%). Table [2] shows 
the average distance between the centres of the real structures 
and the centres of their detected counterparts. The density peaks 
allow to identify the positions of real groups with a good accu- 
racy. 

We also evaluate the ability of the algorithm to separate real 
structures that are very close both in redshift and angular posi- 
tion. In Table [3] we present, for different intracluster distances, 
the density level at which couples of real groups appear as sep- 
arated peaks. Both at low and high redshift it is not possible to 
separate structures whose centres are closer than 1.0 Mpc on the 
plane of the sky and 2<x, in redshift. For larger separations, using 
higher thresholds (5 or 6 cr above the average p) it is possible to 
separate the groups. 

5. A catalogue of the detected overdensities in the 
GOODS-South field 

An inspection of the 3-D density map shows some complex high 
density structures distributed over the entire GOODS field. In 
particular, we find diffuse overdensities at z ~ 0.7, at z ~ 1, 
at z ~ 1.6 and at z ~ 2.3. S ome of these have already been par- 
tially described in lite r ature dGilli et al.l2003|;lAdami et al.l2 005 ; 
Vanzella et al. 200 51 [Trevese et al] 120071: iDfaz-Sanchez et al] 
2007; Castellano et al. 2007). Fig.Q]shows the position of these 
overdensities over the photometric redshift distributions of our 
sample. These overdensities are also traced by the distribution 
of the spectroscopically confirmed AGNs in our catalogue, as 
shown in the lower panel of Fig.[T](these objects are not included 
in the sample used for the density estimation). This link between 
large scale structures and AG N distribution wa s already noted, 
at lower redshift in the CDFS dGilli et al.ll2003l), in the E-CD FS 
dSilverman et alj|2008l) and in the CDFN dBarger et al]|2003l) . 

Within these large scale overdensities, we identify the struc- 
tures, with the procedure described in Sect. [4] Using an analysis 
with a 5cr threshold, we find that two structures identified with 
p > p + 4cr, at z ~ 0.7 and z ~ 1, are the sum of two different 
structures, so we used a 5cr threshold to separate these peaks. 
We then associate the galaxies belonging to the region of over- 
lap between the two structures to the less distant peak. 

Overall, we find four structures at z ~ 0.7, four structures 
at z ~ 1, one at z ~ 1.6 (see also C07) and three structures at 
z ~ 2.3. The density isosurfaces of the structures at z ~ 0.7, 
at z ~ 0.96, z ~ 1.05 and at z ~ 2.3 are shown in Fig. [2j su- 
perimposed on the ACS zsso band image of the GOODS-South. 
T he analogous image for the overdensity at z ~ 1.6 is showed 
in lCas tellano et al. (2007). In the figure, we indicate with a cross 
the peak position of the identified structures. Other overdensities 
present did not pass our selection criteria described in Sect. [4] 
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Fig. 1. Upper panel: photometric redshift distribution of our 
sample (continuous line). Vertical lines mark the redshifts of the 
detected structures. Lower panel: redshift distribution of spec- 
troscopically selected AGNs in the GOODS-South field (contin- 
uous line); the dashed-line histogram is the distribution of the 
AGNs associated with the overdensity peaks in Table [4] 



All the structures are presented in Table [H where we list the 
following properties: 
Column 1: ID number. 

Column 2-4: The position of the density peak (redshift, RA and 
DEC) obtained with our 3-D photometric analysis. 
Column 5: The number of the objects associated with each struc- 
ture as defined above. This number gives a hint on the richness 
of the structure; however it should not be used to compare struc- 
tures at different redshifts because of the different magnitude in- 
tervals sampled. 

Column 6: The average number of field objects present in a vol- 
ume equal to that associated to the structure, at the relevant red- 
shift. We ca lculated this number by integrating the evolutive LFs 
obtained by Salimbeni et al. (2008). In particular, we integrated 
the LF up to an absolute limiting magnitude calculated using the 
average K- and evolutionary corrections and zsso limiting ob- 
served magnitude as done in Sect. |4] In this way we take into 
account the selection effects given by the magnitude cut in our 
catalogue, as a function of redshift. 

Column 7-8: The M200 and r2oo (assuming bias factors 1 and 2). 
The mass M200 is defined as the mass inside th e radius corre- 
spond ing to a density contrast 5 m = 6 ga i/b ~ 200 (Carlber g et al] 
1997), where b is the bias factor. To estimate the 3D galaxy den- 
sity contrast 6 ga i we count the objects in the photometric red- 
shift range occupied by the structure as a function of the cluster- 
centric radius. We then perform a statistical subtraction of the 
background/foreground field galaxies, using an area at least 2.5 
Mpc (comoving) away from the center of every cluster in the rel- 
evant redshift interval. Finally, the density contrast is computed 
assuming spherical symmetry of the structure. The mass inside 
a volume V of density contrast S ga i is determined adapting to 



Table 4. Overdensities in the GOODS-South field. 



ID 


Redshift 


RA 


DEC 


Members 


Field 


M, 00 b=2-l 


i'20o b=2-l 


Peak Overdensity 


Numb. 


Spectroscopic 




r V ir 


M vir 






(J2000) 


(J2000) 






(M o /10 14 ) 


(Mpc) 


a 


Of Z spe 


position 


(km s^ 1 ) 


(Mpc) 


(M o /10 14 ) 


1 


0.66 


53.1623 


-27.7913 


19 


6 


0.15- 0.3 


0.9-1.1 


7 


6 


0.665 ± 0.001 


446 ± 180 


0.85 


1.0 


2" 


0.66 


53.0630 


-27.8280 


50 


17 


0.2-0.4 


1.1-1.3 


10 












3 


0.69 


53.1690 


-27.8747 


54 


20 


0.3 - 0.5 


1.1-1.4 


6 














0.71 


53.0797 


-27.7920 


92 


37 


0.9 - 3.0 


1.7-2.4 


10 


36 


0.734 ± 0.001 


634 ± 107 


1.3 


3.2 


5 


0.96 


53.0843 


-27.9020 


32 


14 


* 


* 


6 












& 


1.04 


53.0570 


-27.7693 


57 


31 


0.5- 1.1 


1.4-1.8 


8 












7 


1.04 


53.1577 


-27.7660 


60 


26 


0.4 - 0.8 


1.2-1.6 


6 












8 w 


1.06 


53.0697 


-27.8773 


38 


18 


0.2 - 0.5 


1.1-1.3 


10 


6 


1.0974 + 0.0015 


446 ± 143 


0.8 


0.8 


9 e 


1.61 


53.1270 


-27.7140 


50 


24 


2.0-4.9 


2.1-2.9 


7 


6 


1.610 ±0.001 


482 ±217 


1.4 


1.46 


10 


2.23 


53.0763 


-27.7060 


20 


8 


0.8- 1.4 


1.5-1.9 


6 












11 


2.28 


53.1470 


-27.7087 


23 


12 


0.6- 1.3 


1.5-1.8 


10 












12 


2.28 


53.0970 


-27.7640 


19 


7 


0.6- 1.6 


1.3-2.0 


9 













a - Gilli et al. (20031 
b - Adami et al. (2005) 
c - Trevese et al. (2007) 
d - Diaz-Sanchez et al. (2007) 
e - Castellano et al. (2007) 

* We do not present A/200 an d r 2W f° r structure 5 because it is located on the edge of the field and it is very close to structure 8, see fig|2l 



S. Salimbeni et al.: Large Scale Structures in the GOODS-SOUTH Field 



7 




Fig. 2. Density isosurfaces for structures at z ~ 0.7 (a), at z ~ 0.95 (b), z ~ 1.05 (c) and at z ~ 2.3 (d) (average, average +2cr, 
average +3cr to average -HlOcr) superimposed on the ACS zsso band image of the GOODS-South field. Yellow crosses indicate the 
density peak o f each structure, the nu mber is the ID of the structure in Tabl^H For the analogous picture regarding the cluster ID=9 
at z ~ 1.6 see Cast ellano et al.l d2007l) . Other overdensities present did not pass our selection criteria described in Sect. [4] 
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our case the metho d used for spectroscopic data at higher z by 
ISteidel et al.l ( fT998h : 



M = p u -V-(\+5 m ), 



(2) 



in which p u is the average density of the Universe and S m is the 
total mass density contrast related to the galaxy number density 
contrast through a bias factor: 1 + 6 m = 1 + 6 ga i/b. We assu me a 
bias factor b in the range 1 < b < 2 (see Arnout s et al.ll 19991) . 
Column 9: The level of the density peak, measured in number of 
cr above the average volume density. 

We then searched the availab l e spectroscopic publ i c data 
dWolf et all l200lt iLe Fevre et all 12004): ISzokolv et al.l l2004t 
Misn oli et al.ll20Q5l: IVanzella et all 120051 120061 l2008h to check 
if any of the members of the structures have spectroscopic red- 
shifts, to estimate their location and the velocity dispersion when 
the statistics is sufficient. Spectroscopic galaxies are considered 
members of a structure if their redshift was within 4-500km/s 
(i.e. three times the velocity dispersion of a rich cluster) from the 
mode of its redshift distribution. From these data, we estimated 
the average redshift of the structure and the velocity dispersion 
using the biweig ht estimators, computed using the Rostat pack- 
age dBeers et al.lll990l) . with 68% confidence uncertainties ob- 
tained from a Jackknife analysis. 

In Table[5]we present a value for the X-ray count rate in the 
band 0.3-4 kev, the corresponding flux (in the interval 0.5-2 keV) 
and the rest-fra me luminosity ( 0.1-2.4 keV), from the Chandra 
2Ms exposure dLuo et ai]|2008l) . We measure the count rates in 
a square of side of ~ 30", centred on the position of the peak of 
each structure. For the count-rate to fl ux conversion we assume 
as spectrum a Raymond-Smith model dRavmond & Sm ith 1977) 
with T=l keV and 3 keV and metallicity of 0.2 Z . 



Table 5. X-ray observations. 



ID 



Count Rate 
0.3-4 keV 

(10~ 5 ) 



Flux" 
0.5-2 keV 

-16 



0.1 



1 a 

■2.4 keV 



S/N'' 



(10- 16 erg s' 1 cm- 1 ) (10 4J erg s~') 



1 


8.49 


6.80-9.01 


0.12- 0.26 


u.l. 


2 


5.56 


4.45-5.90 


0.08- 0.18 


u.l. 


3 


10.1 


8.15-10.98 


0.16- 0.37 


u.l. 


4 


11.2 


9.04-12.31 


0.19- 0.44 


u.l. 


5 


23.7 


19.31-29.21 


0.86- 2.36 


11.3 


6 


5.90 


3.04-4.14 


0.26- 0.76 


u.l. 


7 


5.77 


2.97-4.05 


0.26- 0.74 


u.l. 


8 


9.88 


5.10-6.91 


0.47- 1.37 


u.l. 


9 


5.68 


3.08-4.14 


0.83- 3.67 


u.l. 


10 


9.37 


5.39-7.54 


3.50- 22.43 


u.l. 


11 


5.72 


3.29-4.69 


2.27- 15.06 


u.l. 


12 


6.70 


3.85-5.50 


2.66- 17.64 


u.l. 



a - Values for a Raymond-Smith model with assumed temperature 
respectevelly of 3 keV and 1 keV and metallicity 0.2 Z Q . 
b - u.l. indicates structures with a 3 cr upper limit in the flux. 



5.1. Structures atz~ 0.7 

At redshift z ~ 0.67 we isolate three high density peaks (ID=1,2 
and 3) tha t are part of a larg e scale structure already noted, as a 



ana J ) ma t are part o t a large 
whole, bv lGilii et al.l (l2003h . 



For the structure with ID=1, we estimate the redshift from 
the available 6 spectroscopic data. We find an average red- 
shift of 0.665 ± 0.001 and a velocity dispersion of 446 + 
180 km s" 1 . Assuming that the cluster is virialised, we estimate 
r vir = 0.8Mpc and M vir =1.0- 1O 14 M , using the relations in 
iGirardi et al.ldl998bl) . This estimate is based also on the assump- 
tion t hat there are no infa lling galaxies and that the surface term 
(e. g. ICarlberg et alJll996T) is negligible. Considering the uncer- 
tainties, also due to the small number of spectroscopic galaxies, 
M V! > is fairly consistent with the M200 estimated from the galaxy 
density contrast (0.9 - 3 • 10 14 M o ). 

We also derive the upper limits on the X-ray luminosities for 
this structure, that is of the order of 0.2 - 0.3 ■ 10 43 erg s . All 
the properties presented are consistent wit h the structure being a 
galaxy group/small cluster (Bahcall 1993). 

The structures with ID=2, 3 have upper limits on their X- 
ray luminosities of the order of 0.2 - 0.3 • 10 43 erg s~ l , and their 
masses are of the order of M200 ~ 0.2-0.5- 10 14 M G . These X-ray 
luminositi es and mass es are all typical of galaxy groups/small 
clusters dBahcallll999l) . Each of these structures contains a spec- 
troscop ically confirmed g alaxy detected in the VLA 1.4 GHz 
survey (Mill er et alJl2008n . 

At a slightly higher redshift (z ~ 0.7) we identify a high 
density peak (ID = 4) embedded in another large scale struc- 
ture which was already known in lite rature dGilli et al.| [2003: 
Ada mi et all 120051: iTrevese et all 120071) . In our previous paper 
dTrevese et alj [2007) we identify this structure applying our al- 
gorithm to the data from the K20 catalogue, and classified it as 
an Abell cluster. 

In this new analysis we find that this structure is symmet- 
ric and has a regular mass profile. It has 92 associated ob- 
jects (M B (AB) < -18), and two AGNs. From the density con- 
trast we obtain an r2oo = 1.7 - 2 AM pc and a total mass of 
M 200 = 0.9 - 3.0 ■ 1O 14 M for bias factor b=2-l. From the 
36 galaxies with spectroscopic redshifts, we estimate a red- 
shift location of 0.734 + 0.001 and a velocity dispersion of 
634 + 107km s -1 . We derive a virial radius r vir = \3Mpc, and a 
virial mass M„> = 3.2 ■ 10 14 M G , in good agreement with M2oo- 
The 3 sigma upper limit for the X-ray luminosity in the inter- 
val 0.1-2.4 keV is very low (L x = 0.19 - 0.44 • 10 43 erg s' 1 ). 
Note that the area we considere d does not include the X-ray 
source 173 of iLuo et alj d2008l) . that, similarly to iGilli et al.l 
(2003), we associate to the halo of the brightest clu s ter ga laxy 
(1£>goods-music=9792). Alternatively. lAdami et all d2005l) as- 
sociated the bolometric luminosity (Lx = 0.11 • 10 43 erg s ) of 
the X-ray source 173 to the thermal emission of the intra-cluster 
medium (ICM). From this value they deduced a galaxy velocity 
dispersion around 200 - 300 km s -1 . This value is apparently in 
contrast with the cr,, estimated from the spectroscopic redshifts. 
We also associate to the galaxy IDcoods -music =9792 the object 
236 detected in the VLA 1 .4 GHz survey. It ha s an integrated 
emission of 517.5 + 13.1^ Jv dMiller et alj |2008h. 

From this analysis we can conclude that our two independent 
mass estimates (M200 and M vl >) are consistent with this struc- 
ture being a virialised poor cluster. However, the X-ray emis- 
sion is significantly lower than what is expected from its optical 
properties, as it shows from the comparison in Fig. [3] with the 
M200-LX relations f ound bv lReiprich & Bohringerl d2002l) and by 
Rvko ffetalJd2008l) . 



5.2. Structures at z ~ 1 

At redshift ~ 1 we find four structures (ID= 5, 6, 7 and 8). 
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Fig. 3. Lx vs M200 for the clusters in Tab [5] The horizontal error 
bar is calculated considering a bias factor in the range 1 < b < 2, 
while the vertical error bars are computed varying the gas tem- 
perature between T=l keV and T=3 keV as discussed in the 
main body. The clusters at z ~ 0.7 and z ~ 1.6 are indicated 
by red points and error bars. The M 200 - Lx relat i ons fo und by 
iReiprich & Bohringerl (|2002|) and bv lRvkoff etail (120081) are in- 
dicated by a black and green line respectively. 



The structure with ID=5 at z ~ 0.96 has 32 member galaxies. 
This structure can be associate d to the ex t ended X-ray source 
number 183 in the catalogue by iLuo et alj (120081) derived from 
the 2MS Chandra observation. This extended X-ray source had 
not been associated to any structure so far. From the count rate in 
the interval 0.3-4 keV (S/N=l 1.3) we estimate a luminosity Lx = 
0.86 - 2.36 • 10 43 erg s~ l (in the interval 0.1-2.4 keV). For the 
structures with ID=6, 7 we estimate rioo ~ 1.2— \&Mpc, and a 
total mass ofM2oo = 0.4- 1.1 • 1O 14 M . The 3 sigma upper limits 
for their X-ray luminosity are all slightly below 10 43 erg s~ l , 
consistent with their M200 masses. 

The structure with ID=8 atz ~ 1.06, has 38 associated galax- 
ies, and an AGN spectroscopically confirmed. We derive a pre- 
cise redshift location of z — 1 .0974 + 0.0015 and a velocity dis- 
persion of 446 + I43kmsec~ 1 , from 6 galaxies with spectroscopic 
redshift. From these galaxies we also obtain M vir = 0.8 • 1O 14 M 
and r V i r = O.SMpc. We estimate rioo = 1.1 — \3Mpc, and 
Miqq = 0.2 - 0.5 ■ 1O 14 M , which are compatible values with 
a group of such M„> and ry ir . This structure w as already found 
with different methods by Adam i~et~aLl (120051) . using a friend- 
of-friend algorithm on spectroscopic data fr om the VIMOS VLT 
survey (structure 15 in their Table 4), and bv lDfaz-Sanchez et al.l 
d2007l) studying the extremely red objects on GOODS-South 
(they call this structure GCL J0332. 2-2752). Their redshift posi- 
tions and the velocity dispersions are consistent with those ob- 
tained in the present analysis. The 3 sigma upper limit for the 
X-ray luminosity is around 10 43 erg consistently with the 
estimated M200 mass. 



Considering their properties, these four structures can be 
classified as groups of galaxies. Consistent results fo r the struc- 
ture with ID =6 were obtained in lTrevese et al.l (120071) . 

5.3. Structures at high z 

At redshift z ~ 1.6, we find a compact structure that corre- 
sponds to a f orming cluster, a s already discussed in detail by 
C07 (see also [Kurk et al 2008). We find a regular mass profile 
for this structure, and we estimate an ^200 = 2.1 - 2.9Mpc, and 
a M 2 oo = 2.0 - 4.9 • 1O 14 M . This structure has 50 members, in- 
cluding 3 spectroscopic redshifts, and a confirmed AGN, from 
the GOODS-MUSIC catalogue. We add three other spectro - 
scopic redshifts from the GMASS sample (ICimatti et al.l l2008). 
From these 6 redshifts we estimated a velocity dispersion of 
482 + Hlkmls, and derived an M vir = 1.4 • 1O 14 M and 
r V i r - \A6Mpc. This estimate is consistent with the value in 
Table [4] We derive an upper limit to the X-ray luminosity of 
0.83 - 3.67 ■ 10 43 erg s' 1 (0.1-2.4 KeV), lower than expected 
from the velocity dispersion and the estimated M200 (see Fig . [3j . 

At z ~ 2.2 we find a diffuse overdensity, similar to those 
at lower redshift, embedding three structures. We associate to 
these structures 20, 23 and 19 galaxies. We estimate for all these 
structure an rim ~ 1.3 - 2Mpc and a mass of M200 ~ 0.6 - 1.6 • 
10 14 M . These structures appear to be comparable to those at 
~ 0.7 and ~ 1.6, and they could be forming clusters. 



6. Colour-Magnitude diagrams 

We study the colour magnitude diagrams (U - B vs. Mb) for all 
the structures, as shown in Figs. [4] and [5] To estimate the slope 
of the red-sequence, we define its members as passively evolving 
galaxies according to the physical criterion age It > 4, where the 
age and t (the star formation e-folding time) are inferred for each 
galaxy from the SED fitting (Sect. |2). This quantity is, in prac- 
tice, the inverse of the Scalo parameter (lScalolll9 86) and a ratio 
of 4 is chosen to distinguish galaxies having prevalently evolved 
stellar populations from galaxies with recent episodes of star for- 
mation. Indeed, an age It = 4 corresponds to a residual 2% of the 
initial SFR, f or an exponent i al star formation history, as adopted 
in this paper. iGrazian et al.1 d2006bl) showed that this value can 
be used to effectively separate star forming galaxies fro m the 
passively evolving population (see IGrazian et al.l l2006bl also 
for the discussion on the uncertainty associated to this parame- 
ter). Passively evolving galaxies are indicated in figures as filled 
squares. 

Fig. [4] shows the colour magnitude diagrams for the four 
structures between z — 0.66 andz = 0.71. The cluster at z ~ 0.71 
(Panel d) shows a well defined red sequence, while the three 
structures at z ~ 0.66 have fewer passively evolving galaxies. 
Therefore, in order to increase our statistics, we estimate the 
colour-magnitude slope combining all the four structures in the 
interval 0.66 < z < 0.71 (see Panel a in Fig . [5j ■ We obtain a value 
-0.023 + 0.006 for the slope. The resulting colour-magnitude re- 
lation is plotted in all panels in Fig. [4] and in the panel a in Fig. 
[5] as a continuous line. The dotted lines constrain the error at 1- 
sigma obtained with a Jackknife analysis. It is possible to see 
in Fig.|4]that this average colour magnitude relation is roughly 
consistent with the position in the (U-B) vs B diagram of the 
galaxies belonging to each single structure. We therefore apply 
the same method at higher redshift, i. e. we estimate the slope 
of the red sequence by combining the different structures at the 
same redshift. 
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Fig. 4. Rest frame colour magnitude relations (U - B vs M^) for 
each structure at z ~ 0.7. Squares indicate passively evolving 
galaxies selected as age It > 4, and the circles are galaxies with 
age/r < 4. Filled points indicate galaxies with spectroscopic 
redshift. The continuous lines are the fit to the red sequence of 
all the combined structures. The dotted lines are the uncertainties 
at 1-sigma obtained with a Jackknife analysis. 



Fig- El panel b, shows the colour magnitude diagram for the 
structures at z ~ 1. We find a slope of -0.03 + 0.01. 

Panel c in Fig. [5] shows the colour magnitude diagram for 
the structure at z ~ 1.6. In this case we have galaxies distributed 
on less than a magnitude range, that is insufficient to estimate 
the slope of the "red sequence". However, if we plot the two 
sequences obtained at lower redshift, we can see that the few 
passively evolving galaxies are consistent with them. 

Finally, at redshift ~ 2, we have only 4 passive objects from 
the combination of 3 structures and there is no evidence of a well 
defined red sequence. We note that the colours of these objects 
are generally bluer in comparison to the colour of the relations 
found at lower redshifts. 

The values of the slopes of the structures at redshift ~ 0.7 
and ~ 1 are consistent wi t h those of previous d eterminations 
(e.gJBlakeslee et alj 120031; iHomeier etaf] 120061: pTrevese et alj 
2007). We confirm that the observations indicate no evolution 
up to redshift ~ 1. This would imply t hat the mass-metal licity 
relation that produces the red sequence ( Kodam aet alJI 1998b re- 
mains practically constant up to, at least, z ~ 1. 



Fig. 5. Same as Fig. [4] Panel a: structures at z ~ 0.7; panel b: 
structures at z ~ 1; panel c: structure at z ~ 1.6 ; panel d: struc- 
tures at z ~ 2.2. The two dashed lines in each of the two last bins 
of redshift are the red sequences estimated at z ~ 0.7 (shallower 
slope) and ~ 1 (steeper slope) 

7. Galaxy properties as a function of the 
environment 

To each object in the sample we associate the comoving density 
at its position, and we study galaxy properties as a continuous 
function of the environmental density. 

7.1. Galaxy populations: bimodality 

We study the variation of the fraction of red and blue galaxies 
as a function of the environmental density. To separate red and 
blue galaxies we use the minimum in the bimodal galaxy distri- 
bution in the (U - V) vs . B colour magnitude diagram, derived by 
Salimbeni et al.ld2008l) . Fig.|6]shows the fraction of red and blue 
galaxies for different rest frame B magnitudes in four redshift in- 
tervals. In general, for every environment, we find that, at fixed 
luminosity, the red fraction increases with decreasing redshift, 
and, at fixed redshift, it increases at increasing B luminosity. We 
also find that for z < 1 .2 the red fraction increases with density 
for every luminosity, while this effect is absent at higher redshift. 

Our result s exte nd to higher redshift those obtained by 
ICucciati et alj (120061) on the VVDS survey, with a shallower 
spectroscopic sample that reaches z ~ 1 .5. We find that at z > 1 .2 
even the highest luminosity gal axies are blue, star fo rming ob- 
jects, similarly to the results in Cucc iati et al.l (|2006), although 
our colour selection is slightly different, since we select in colour 
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Fig. 6. Fraction of red (filled circles) and blue galaxies (filled triangles) at decreasing rest frame B magnitudes (from top to bottom) 
in four contiguous intervals of increasing redshift (from left to right). Vertical errorbars indicate the poissonian uncertainty in each 
bin. The shaded areas are obtained by smoothing the red (blue) fraction with an adaptive sliding box. The horizontal errorbars 
indicate the range of density covered by the 5-95 % of the total sample. 



two complementary samples, while they select two extreme red 
and blue populations {{u* - g') > 1.1 and (u* - g') < 0.55)). Our 
results are also in agreem ent with the analysis of the DEEP2 sur- 
vey bv lCooper et al.l d2007l) in the redshift range 0.4 < z < 1 .35. 
They found a weak correlation between red fraction and density 
at z ~ 1.2. We see that at z > 1.2 this correlation disappears, 
indicating that the change probably occur in the critical range 
1 .5 < z < 2.0, at least in the environments probed by our sample. 
However we note that, given the relatively small area covered, 
we do not probe very high density regions (i.e. rich clusters), at 
variance with wi de, low redshift surv eys. When rich clusters are 
considered (e.g. iBalogh et al.ll2004l) . a stronger variation with 
environment in the colours of faint galaxies is seen. In any case, 
the disappearance at z > 1 .2 of the variation of the red fraction 
in the density range probed by our sample, is an indication that 
a relevant change in galaxy properties takes place at z ~ 1 .5 - 2. 

7.2. Galaxy physical properties in high and low density 
environments 

We then study the distribution of physical parameters and photo- 
metric properties for galaxies in high density environments, and 



compare it to field galaxies. The first sample is defined as the 
combination of the data from structures with similar redshifts 
('group galaxies' hereafter). The field galaxies are defined as 
those with an associated p lower than the median density (0.0126 
for z < 1.8 and 0.0085 for z > 1.8) of the entire sample ('field 
galaxies' hereafter). We quantify the differences in the distribu- 
tions of the galaxy physical properties, i.e. mass, age, star for- 
mation rate, through the probability P KS of the two samples, ob- 
tained as described above, using a Kolmgorov-Smirnov test. We 
reject the hypothesis that two samples are drawn from the same 
distribution if P KS < 5 • 10~ 2 . 

Fig. [7] shows the distribution of the galaxy total stellar mass 
in high and low density regions, in the same four contiguous 
redshift intervals used before. The galaxies in high density envi- 
ronment have a distribution that generally peaks at higher masses 
with respect to "field" galaxies. For the mass distribution we find 
a significant difference in all but the last redshift bin as shown 
from the Pks ■ It is important to remark here that the shape of the 
distributions at low masses could depend on the luminosity se- 
lection. In fact, a magnitude-limited sample does not have a well 
defined limit in stellar mass. This effect depends on the range 
of M/L ratio spanned by galaxies with different colours, e. g. as 
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shown in Font ana et al.l {2006) in our sample, at z ~ 1 , Mj Lk ex- 
tends from 0.9, for redder objects, to 0.046, for bluer objects. If 
a colour segregation is present as a function of the environment, 
it could bias the distribution favouring the observation of lower 
mass galaxies in less dense regions, where the fraction of blue 
galaxies is higher. Although, as shown in Fig. [6] we do not find a 
strong colour segregation, especially at z > 1, we carry out here 
also a more conservative analysis. We consider only the range 
of masses above the completeness mass limit obtained from the 
maximal M/L- S50 for a passive evolving system (log(M) > 9.0 at 
z ~ 0.6, log(M) > 9.6 at z ~ 1, log(M) > 10.5 at z ~ 1.6 and 
log(M) > 11.1 at z ~ 2.15). Considering galaxies above these 
mass limits we find that the masses of "group" galaxies are still 
higher than those of "field" galaxies, for the lower bin in redshift 
(Pks — 9.7 ■ 10~ 4 ). At z > 1 .2, however, it is not possible to give 
a conclusive result due to the low statistic caused by this mass 
cut. 




log(M) log(M) 

Fig. 7. Galaxy stellar mass distribution in four redshift intervals. 
Shaded red histograms represent galaxies associated with the 
density peaks and empty black histograms represent galaxies in 
the low density regions, as described in the text. In each panel the 
average value of log(M) for the two distributions are indicated 
by arrows of the same colour. The K-S probability is reported in 
each panel. Vertical lines indicate the mass limit at the median 
redshift of the bin (log(M) = 9.0 at z ~ 0.6, log(M) = 9.6 at 
Z~ 1, log(M) = 10.5 atz ~ 1.6 andlog(M) = 11.1 atz ~ 2.15). 

Analogous results are found from the analysis of the lumi- 
nosity distribution of "field" and "group" galaxies. In particular, 
we find that the distribution of galaxies in "groups" have on av- 
erage brighter Mj rest-frame magnitudes at all redshifts. The re- 
sults are similar also for the other rest frame bands, implying that 
galaxies in high density environments have, on average, greater 
bolometric luminosity with respect to field galaxies. 

Finally, we study the age and SFR distributions for "group" 
and "field" galaxies (see Fig. |S). Only at low redshift there ap- 
pears to be a significant difference (respectively Pks =3.0- 10~ 2 
and Pks = 6.7 • 10~ 3 , see Fig. [HJ. The two age distributions show 



a similar shape for young galaxies, but "group" galaxies have a 
higher fraction of old galaxies. As also shown by the difference 
in the average ages for the two samples, the "group" galaxies 
are older than the field ones. At higher redshifts the two distri- 
butions do not show significative differences. Indeed, at higher 
redshifts, any possible difference in the age of the two galaxy 
populations is probably smaller than the uncertainty on the ages. 
Analogously, star forming galaxies have a similar distribution 
for "group" and "field" samples, but the "group" sample has a 
higher faction of galaxies with low star formation as it is also 
shown by the different values of the average SFR. 




2 4 6 

Age [Gyr] 




Log(SFR [M -yr->] ) 

Fig. 8. top: As in Fig.|7]but for the ages of galaxies. The average 
value of the age for the two distributions are indicated by arrows. 
Bottom: As in Fig.|7]but for SFR of galaxies. The average values 
of the log(5 FR) for the two distributions are indicated by arrows. 



8. Summary and Conclusions 

We applied a (2+l)D algorithm on the GOODS-MUSIC cata- 
logue to identify structures in this area. This algorithm com- 
bines galaxy angular positions and precise photometric redshifts 
to give an adaptive estimate of the 3D density field effectively 
also at z > 1, and in a wide area. In this way we obtained a den- 
sity map from redshift 0.4 up to 2.5, and we isolated the higher 
density regions. To identify density peaks we chose a conserva- 
tive selection criteria (at least five galaxies in connected regions 
of p > p + Act) in order to maximise the purity of our sample. 

We built mock catalogues simulating the GOODS-South 
field. Applying our density thresholds and selection criteria on 
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these catalogues, we found a purity near 100% (with less than 
15-20% of lost structures) up to redshift 1.8, and ~ 75-80% at 
higher redshift. In the higher redshift range, the criterium is very 
conservative, to keep a low number of false detections, therefore 
the completeness is low (< 40%). From the simulations we also 
evaluated the ability of the algorithm in separating real structures 
that are very close both in redshift and angular position. Both 
at low and high redshift it is not possible to separate structures 
whose centres are closer than 1 .OMpc on the plane of the sky 
and 2cr z in redshift. For larger separations it is possible to distin- 
guish the groups, but using higher thresholds (5 or 6 cr above the 
average). 

We found large scale overdensities at different redshifts (~ 
0.6, ~ 1, ~ 1.61 and ~ 2.2 ), which are well traced by the AGN 
distribution, suggesting that the environment on large scales 
(~ lO Mpc) has an influence on AGN evolution (Silverman et al.l 
2008). We isolated several groups and small clusters embedded 
in these large scale structures. Most of the structures at z ~ 0.7 
and ~ 1 have properties of groups of galaxies: their masses are 
of the order of M 20 o = 0.2 - 0.8 ■ 10 14 M o , and their X-ray lu- 
minosities are slightly below 10 43 erg s~', consistent with the 
expectations of the M2oo-La- relations. The structure at z — 0.71, 
and those at z > 1 .6 seem to be more massive, and in particular 
the structures with ID=4, and 9 can be classified as poor clusters. 
It is interesting to note that both these structures are significantly 
X-ray underluminous, as it is evident by a comparison w ith the 
M2oq-Lx relations f ound bv lReiprich & Bohringerl (l2002l) and by 
Rvkoff et al. (2008) (Fig. [5j. This is not surprising since several 
authors have observed that optically selected structures have an 
X-ray emission lower than what is expected from the observa- 
tions of X-ray selected groups and clusters: th is effect has been 
obser ved at low redshift bot h in small groups ( Rasmusse n~et~aT1 
2006) and in Abel l clusters (Popess o et al . 2007) and in clusters 
at 0.6 < z < 1.1 dLubin et al.ll2004l) . These results may be ex- 
plained by the fact that such optically selected structures are still 
in the process of formation or the result of the alignment of two 
substructures along the line of sight, although it cannot be ex- 
cluded that they contain less intracluster gas t han expected, be- 
cause of the effect of strong galactic feedback (Rasmusse n et al.l 
2006). If these structures are virialised, as probable in the case 
of the massive structure at z = 0.71 (ID=4), this may be an in- 
dication that they contain less intracluster gas than expected. It 
is worth investigating this issue in future deep surveys, since it 
would have interesting implications on the evolution of the bary- 
onic content of these structures. 

We then studied the colour magnitude diagrams (U - B vs 
Mb ) for all the structures. We defined the members of the red- 
sequence according to the physical criterion age It > 4 which 
should select passively evolving galaxies with little residual star 
formation. We confirmed no evolution of the red sequence slope 
up to redshift ~ 1 . This implies that the mass-metallicity relation 
that produces the slope of the red sequence remains constant up 
toz ~ 1. 

We then studied the variation of the fraction of red and blue 
galaxies as a function of the environmental density. We found 
that, at fixed redshift, the red fraction increases at increasing B 
luminosity, while, at fixed luminosity, it increases with decreas- 
ing redshift. We found that the increment of the red fraction at 
growing density disappears at z > 1.2. 

We also studied galaxy properties in different environments. 
We found that the galaxies in high density environments have 
higher masses with respect to "field galaxies", in qualitatively 
agreement with a downsizing scenario. The mass distributions 
show a significant difference in all but the last redshift bin. 



Similarly, the galaxies in groups have on average brighter rest- 
frame magnitudes and there is a greater number of bright galax- 
ies in groups at all redshifts compared to field galaxies. Finally, 
the age and SFR distributions for the two subsamples appear dif- 
ferent only at low redshifts where "group galaxies" are generally 
older and less star forming than"field" ones. 

From the analysis of the environmental dependence of 
galaxy colours and mass as a function of redshift, and from the 
absence of any well defined red sequence at high redshift, we 
can argue that a critical period in which some basic characteris- 
tics of galaxy populations are established is that between z ~ 1.5 
and z ~ 2. 
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